UNCLASSIFIED 


_ AD  NUMBER _ 

AD879556 

LIMITATION  CHANGES 
TO: 

Approved  for  public  release;  distribution  is 
unlimited. 


FROM: 

Distribution  authorized  to  U.S.  Gov't,  agencies 
only;  Test  and  Evaluation;  JAN  1971.  Other 
requests  shall  be  referred  to  Office  of  Naval 
Research,  875  North  Randolph  Street,  Arlington, 
VA  22203-1995. 


_ AUTHORITY 

ONR  ltr,  Aug  1973 


THIS  PAGE  IS  UNCLASSIFIED 


OKI  *>8795  56 


ADL- 725 80-4 


TECHNICAL  REPORT  4 


PASSIVE  RANGING  AND  TARGET  MOTION  ANALYSIS 


James  M.  Dobbie 


Under  Contract 


N0001 4-70-C-0322 
NR  364-025/2-12-70  (462) 


Prepared  For 


Office  of  Naval  Research 
Department  of  the  Navy 
Washington,  D.  C.  20360 


Prepared  By 


Arthur  D.  Little,  Inc. 
Acorn  Park 

Cambridge,  Massachusetts  02140 


January  1971 


"REPRODUCTION  IN  WHOLE  OR  IN  PART  IS  PERMITTED  FOR 
ANY  PURPOSE  OF  THE  UNITED  STATES  GOVERNMENT." 


D  D  C 

lEJinQETi 


G^ 


i> 


- . .  -  •  •[  mu  illiiteiitiiil  iWrdf 


I 

_ UNCLASSIFIED _ 

DOCUMENT  CONTROL  DATA  -RAD 

^^^^^JStcurlt)^cloiilltcBtloMMItlo^»o4^^kalro*^on«nn4oah£annotollotM»uilib^rntoro4jmho^h^rorol^ogotM^clo*olllod^ 

1.  orioinatino  activity  (Corpora*  author)  la.  uroar  security  classipic ation 

J  Arthur  D.  Little,  Inc.  UNCLASSIFIED 

Acorn  Park  *6.  sroup 

Cambridge,  Massachusetts  02140 

I.  REPORT  T I TLS 


Passive  Ranging  and  Target  Motion  Analysis 

4.  descriptive  NOTES  (Typ*  of  ropott  ond  Inelutlvt  Salts) 

Technical  Report  4 

»•  AUTHOR!*)  (hr* l  namo,  ai 733a  Initial,  TomTnomo) 

James  M.  Dobbie 


a.  MPORT  DATE 

January  1971 

Aa.  CONTRACT  ON  CHANT  NO. 

N00014-70-C-0322 

b.  PROJECT  NO. 


7k.  NO.  OP  REPS 

5 


ADL- 72580-4 


fsk.  OTf 
•hit 


REPORT  NOIS)  (Any  al har  nuaikara  Mat  atap  fca  aaaljna* 

Technical  Report  4 


to.  DISTRIBUTION  STATEMENT 

Each  transmittal  of  this  document  outside 
must  have  prior  approval  of  the  Office  of 
Underwater  Systems  Center/Newport 

the  agencies  of  the  U.S.  Government 

Naval  Research  Code  462  or  the  Naval 

It.  SUPPLEMENTARY  NOTES 

tl.  SPONSORINS  MILITARY  ACTIVITY 

Naval  Analysis  Programs 

Office  of  Naval  Research 

Dept,  of  Navy,  Washington  D.C.  20360 

"'■  A  review  is  made  of  various  methods  for  passive  ranging  and  target  motion  analysis, 
including  some  new  methods  based  on  three  tracking  legs.  Estimators  for  range  and 
other  parameters  are  compared  with  respect  to  their  responses  to  errors  from  four 
major  sources.  A  program  for  additional  tests  and  comparisons  is  outlined.  /  . 


DD  ,'r..1473 


REPLACES  DO  PORM  1*7*.  I  JAN  14.  WMICN  IS 
OBSOLETE  POR  ARMY  USB. 


UNCLASSIFIED 

Security  Classification 


f 


1 


PASSIVE  RANGING  AND  TARGET  MOTION  ANALYSIS* 

Introduction 

Passive  ranging  and  target  motion  analysis  have  been  studied  for 
many  years.  A  few  of  the  many  methods  that  have  been  developed 
have  been  adopted  for  use.  The  methods  in  use  have  not  been 
satisfactory  in  all  respects,  and  the  search  continues. 

There  appears  to  be  no  "best"  way  to  proceed.  That  is,  there  is 
no  method  that  is  superior  to  all  other  methods  in  all  respects. 
Under  these  conditions  the  advocates  of  various  methods  argue 
their  cases  by  the  adversary  procedure.  A  better  procedure 
would  be  to  examine  the  major  characteristics  of  the  leading 
contenders  on  an  objective  basis,  with  the  same  tests  and  condi¬ 
tions  applied  to  all.  The  main  question  remaining  then  would  be 
that  of  the  characteristics  that  are  important  in  a  particular 
application,  and  the  relative  weights  that  should  be  assigned  to 
the  characteristics  in  choosing  an  estimator  for  the  range  or 
for  a  complete  set  of  target  motion  parameters. 

The  purpose  of  this  report  is  to  describe  some  of  the  leading 
contenders,  to  discuss  their  characteristics  to  the  extent  that 
they  are  known,  and  to  outline  a  test  and  comparison  program. 

Approach 

The  major  characteristics  of  a  range  estimator,  or  an  estimator 
of  some  other  parameter,  are  the  following: 

(a)  Formula  accuracy 

(b)  Sensitivity  to  error  measurements 

(c)  Ruggedness  to  target  maneuvers 


This  research  was  supported  by  the  Office  of  Naval  Research 
under  Contract  No.  N00014-70-C-P322 . 
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All  three  characteristics  are  important  in  most  applications, 
but  the  relative  weights  may  be  different. 

By  formula  accuracy  we  mean  the  accuracy  of  the  estimator  when 
there  are  no  measurement  errors  of  any  kind  and  the  target  motion 
is  strictly  linear.  The  reader  who  has  a  classroom  knowledge  of 
statistics  and  little  experience  with  estimation  problems  of  the 
type  we  are  considering  may  be  surprised  that  we  would  consider 
any  estimators  that  do  not  have  formula  accuracy.  Yet,  as  will 
be  seen  later,  some  estimators  of  this  type  are  not  only  useful 
but  may  be  the  best  we  can  do  in  regard  to  certain  character¬ 
istics. 

The  above  remarks  are  made  to  dispel  any  doubt  that  the  estima¬ 
tion  of  range  and  range  rate  from  bearings  only  is  a  difficult 
problem.  The  theory  of  unbiased  estimators  having  minimum 
variance  just  doesn't  apply.  All  estimators  of  range  and  range 
rate  that  we  know  of  are  biased.  Some  are  biased  through  for¬ 
mula  inaccuracies,  and  all  are  biased  stochastically. 

Almost  all  estimators  are  derived  under  the  assumption  of  linear 
target  motion.  This  assumption  in  the  derivation  of  possible 
estimators  is  justified,  since  bearing  accuracies  are  not  suf¬ 
ficient  to  predict  non-linear  target  motion.  That  is,  the  best 
that  we  can  do  with  present  accuracies  is  to  predict  the  linear 
component  of  target  motion,  and  even  this  much  is  in  doubt,  in 
my  opinion.  But,  whatever  assumptions  are  made  in  the  derivation, 
we  must  recognize  the  fact  that  the  target  may,  and  often  does, 
have  motion  that  is  far  from  linear.  Hence,  the  estimator, 
whether  derived  from  the  linear  assumption  or  not,  should  be 
subjected  to  tests  against  non-linear  motion. 

Wien  we  consider  the  difficulties  of  the  problem  and  the  several 
likely  sources  of  large  biases,  it  seems  evident  that  highly 
sophisticated  techniques  are  not  justified.  At  least  they  should 
not  be  used  until  it  can  be  demonstrated  that  they  significantly 
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reduce  the  total  bias.  The  elimination  of  a  random  component  of 
the  error  in  the  range  estimate,  for  example,  will  do  little  to 
increase  the  effectiveness  of  the  estimate,  if  the  standard 
deviation  of  the  random  component  is  small  relative  to  the  total 
bias.  A  significant  improvement  in  the  estimate  requires  a  sig¬ 
nificant  reduction  in  the  dominant  error,  in  general.  And,  the 
dominant  error  usually  is  a  bias. 

There  are  four  sources  of  bias,  as  follows: 

(1)  Inaccuracies  in  the  formulas 

(2)  Biases  in  the  bearing  observations,  or  other 
observations 

(3)  Stochastic  bias 

(4)  Non-linear  target  motion 

The  first  one  has  been  discussed  above.  Biases  in  the  bearing 
observations  are  known  to  exist.  The  bias  depends  on  the  angle 
off  the  bow  of  the  tracking  submarine,  and  usually  is  called  the 
delta  bias.  The  third  and  fourth  sources  of  bias  also  have  been 
mentioned  above.  We  believe  that  all  four  sources  of  bias  must 
be  included  in  an  adequate  analysis. 

If  the  target  maneuvers  radically  during  the  tracking  interval, 
the  fourth  bias  almost  certainly  is  the  dominant  bias.  If  the 
maneuver  is  mild,  one  of  the  other  biases  may  be  the  dominant 
one. 

The  net  effect  of  the  biases  and  random  errors  will  depend  on 
the  use  that  is  made  of  the  estimates,  and  the  measure  of  per¬ 
formance  that  is  chosen.  In  general,  a  reduction  in  the  magni¬ 
tude  of  the  dominant  bias  component  will  produce  a  net  gain, 
while  a  reduction  in  the  dominant  variance  may  produce  a  net 
gain  or  a  net  loss.  To  avoid  considering  various  applications 
and  measures  of  performance  we  will  concentrate  on  the  bias 
components,  for  which  the  favorable  direction  is  known. 

3 


Arthur  D  Little  Inc 


t 


For  some  purposes  a  range  estimate  is  sufficient.  For  other 
purposes  a  complete  TMA  is  required.  The  parameter  estimators 
in  a  complete  TMA  may  include  a  range  estimator  directly,  or 
require  an  additional  calculation.  Complete  TMA  estimators  and 
range-only  estimators  are  included  in  this  preliminary  compari¬ 
son  of  bias  components. 


Notation 

We  use  the  standard  naval  coordinate  system.  For  the  ith  obser¬ 
vation  let 


t^  =  time 

b^  ■  bearing  of  the  target  from  the  tracking  submarine 

r^  *  range  of  the  target  from  the  tracking  submarine 

u  ^  ■  range  component  of  target  velocity  measured  outward 

u^  ■  normal  component  of  target  velocity  measured  clockwise 

v^  **  range  component  of  own-submarine  velocity  measured 
outward 

v.^  *  normal  component  of  own-submarine  velocity  measured 
clockwise 


It  will  be  convenient  to  use 


B^  ■  bearing  observation,  or  estimate  of  b^ 


'ij 


"i 

-  t.-t 


estimate  of  r, 


J  i  ’  ij 


VBi  ’  Sij 


sin  B±j  ,  C1J  -  cos  B 


Estimates  From  Two  Tracking  Legs 

We  assume  that  the  tracking  submarine  uses  two  essentially  linear 
tracking  legs  on  two  quite  different  courses,  of  the  type  described 
as  "Lag-lead",  for  example.  Four  bearing  observations  are  suffi¬ 
cient  for  a  complete  TMA,  on  the  assumption  of  linear  target 
motion.  The  corresponding  estimator  for  the  range  r^  at  the 

time  t,  of  the  fourth  observation  is 
4 

(1)  =  t14t24t34(vblS23+Vb2S3l'l'Vb3S12)  ... 

4  t12t34S13S24  "  t13t24S12S34 
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This  equation  is  equivalent  to  equation  (16)  of  reference  (a) . 

It  is  the  only  exact  estimator  for  r^  from  four  bearing  observa¬ 
tions  with  the  linear  assumption.  Estimators  for  other  parameters 
sufficient  for  a  complete  TMA,  also  can  be  written. 


Some  inexact  estimators  for  r^  are  the  following: 
(2) 

R'  "  (vb3  ~  vblC13  ‘  "  334* 


4 

,(3) 


(vb3  _  VblC13)/(B12  ~  634) 


where 


312  “  (S12/t12)(C13/C24)*  3 34  “  S34^34 


(2) 

(3) 

(4) 


Equation  (2)  is  equivalent  to  equation  (6)  of  reference  (a) ,  and 
is  obtained  by  assuming  that  the  u^  component  of  target  motion 
during  the  interval  (t^t^)  is  negligible.  Equation  (3)  is 
equivalent  to  equation  (2-23)  of  reference  (b),  and  is  obtained 
from  equation  (2)  by  assuming  that  the  v^  component  of  own-ship 
motion  during  (t2,t^)  also  is  negligible. 


A  form  of  the  Ekelund  range  estimator  is  obtained  from  equation 
(3)  by  letting  t2  approach  t^  t^  approach  t3>  and  approximating 
and  by  1.0. 


The  estimator  is 

,(*) 


R 


(vb3  ■  Vbl)/(B1  -  V 


(5) 


•  • 

where  and  are  estimates  of  the  bearing  rates  at  t^  and  t^. 

This  estimator  usually  is  a  better  estimator  for  r„  than  for  r, . 

3  4 

A  slightly  better  estimator  for  is 


R 


(5) 


<vb4  -  Vb2)/(B2  ’  V 


(6) 


which  can  be  derived  directly  from  the  assumption  that  the  u 


and  v^  components  are  negligible;  see  reference  (c). 
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Other  forms  similar  to  those  in  equations  (2),  (3),  (5)  and  (6) 
also  have  been  used.  In  particular,  the  special  form  of  equation 
(3)  when  t^  *  t^  is  called  the  DEVGROUP  range  estimate  in  refer¬ 
ence  (b) . 


A  complete  TMA  can  be  obtained  in  several  ways.  First,  we  can 
use  the  exact  solution  based  on  four  bearings.  The  range  compo¬ 
nent  of  this  solution  is  obtained  from  equation  (1),  and  esti¬ 
mators  for  other  parameters  can  be  written.  However,  we  do  not 
recommend  it.  The  range  estimator  is  very  sensitive  to 

random  bearing  errors  and  to  mild  deviations  of  the  target  from 

a  linear  course;  and  the  estimator  for  u  .  is  even  more  so.  The 

r4 

difficulty  stems  from  the  fact  that  the  denominator  in  (1)  is  a 
difference  of  two  terms,  which,  for  a  standard  approach,  have 
the  same  sign  and  order  of  magnitude.  Often  the  magnitude  of 
the  difference  between  the  terms  is  small  relative  to  the  magni¬ 
tudes  of  the  terms.  Hence,  a  small  fractional  error  in  a  term 
produces  a  large  fractional  error  in  the  difference. 


A  simple  TMA  based  on  four  bearing  estimates  is  the  following: 
TMA  No.  1 

=  estimate  of  bearing  at  t  =  t^ 

R^  *  estimate  of  range  at  t  =  t^,  from  equation  (2),  (3), 
or  (5) 

U  .  =  estimate  of  u  .  =0 
r4  r4 

Ub4  =  estimate  of  u^  =  vfa4  +  R4  S^/t^  (7) 


The  four  estimates  above  give  us  the  (polar-coordinates)  position 
and  velocity  components  at  t  =  t4>  The  use  of  0  for  Ur4  is  con¬ 
sistent  with  the  assumptions  made  in  obtaining  the  range  esti¬ 
mators.  An  exact  estimator  for  U,  .  is 


b4 


=  v, 


b4 


+  R3  S34/t:34 
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* 


Hence,  the  estimators  Rp\  and  R,^\  which  are  better 

4  4  4 

estimators  for  r_  than  for  r. ,  are  preferred  over  the  estimator 

(5)  J  4 

R^  in  the  estimator  (7)  for  U^. 

Accurate  estimates  of  the  bearing  rates  that  appear  in  the  Ekelund 
range  estimator  can't  be  obtained  from  the  four  (B^,  t^)  pairs 
alone.  We  need  frequent  bearing  observations.  Then  the  bearing 
rates  can  be  estimated  from  the  smoothed  bearing-time  plot,  or 
computed  directly  from  the  bearing  observations  by  using  some 
method  of  smoothing.  A  method  using  exponential  smoothing  in 
reference  (c)  was  found  to  be  satisfactory.  If  a  more  accurate 
estimate  of  the  bearing  rate  can  be  obtained  from  the  bearing¬ 
time  plot  at  a  time  that  is  not  the  end  of  a  tracking  leg, 
and  t^  at  which  the  velocities  and  bearing  rates  are  estimated 
in  equation  (6)  can  be  replaced  by  earlier  times.  The  increase 
in  accuracy  of  the  bearing  rates  may  more  than  offset  the  lag 
in  the  range  estimate. 

If  we  take  many  bearing  observations  on  the  two  tracking  legs, 
other  methods  of  target  motion  analysis,  such  as  CHURN,  can  be 
used.  Also,  various  filtering  and  smoothing  techniques  can  be 
applied.  These  methods  and  techniques  will  be  discussed  briefly. 

Discussion  of  Range  Estimators  From  Two  Tracking  Legs 

The  range  estimator  R^  has  no  bias  from  formula  accuracy,  which 
is  the  only  advantage  it  has  over  the  other  range  estimators. 

The  biases  from  formula  inaccuracies  in  the  other  range  esti¬ 
mators  vary  from  trivial  to  serious,  depending  on  many  factors, 
such  as  the  range  at  the  start  of  the  tracking  interval,  the 
duration  of  the  tracking  interval,  and  the  motions  of  the  target 
and  own  submarine  during  tracking.  Biases  exceeding  15  percent 
are  common,  and  biases  exceeding  30  percent  are  not  unusual. 

Some  results  from  computer  simulations  are  given  in  references 
(b)  and  (c). 
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The  bias  in  the  bearing  observations  usually  is  called  the  delta 
bias.  It  depends  on  the  angle  off  the  bow  of  the  tracking  sub¬ 
marine  at  which  the  bearing  is  taken,  and  it  usually  varies  slowly 
with  changes  in  this  angle.  Hence,  the  delta  bias  is  nearly 
constant  for  a  tracking  leg,  but  varies  from  leg  to  leg.  Hence, 
by  taking  and  on  the  first  leg,  and  and  B^  on  the  second 

leg,  the  errors  in  and  the  bearing  rates  can  be  kept 

small.  Thus,  Rp\  R/^,  and  R.^  have  a  distinct  advan- 

(1)  4  4  4  4 

tage  over  R.*1  for  this  bias  component.  A  partial  analysis  for 
(5)  4 

R^  was  made  in  reference  (c) . 


Random  bearing  errors  also  produce  larger  errors  in  than  in 

the  inexact  estimators,  both  in  bias  and  in  variance.  The  biases 
arise  from  the  fact  that  the  range  estimators  are  not  linear  in 
the  bearing  inputs.  This  source  of  bias  in  the  range  estimator 
probably  is  not  the  dominant  bias  for  the  range  estimators  dis¬ 
played  above,  although  we  do  not  know  of  any  results  that  apply 
to  the  question.  (However,  it  may  be  a  serious  bias  in  the 
CHURN  TMA  under  some  circumstances;  see  the  discussion  of  CHURN 
in  a  later  section.) 


If  the  target  departs  mildly  from  a  linear  course,  say  10  degrees 
change  in  course  near  the  middle  of  the  tracking  interval,  the 
bias  in  R^1^  is  larger  than  chat  in  the  other  estimators.  If  the 
target  maneuvers  radically,  say  in  the  middle  half  of  the  track¬ 
ing  interval,  all  range  estimators  are  useless.  To  avoid  accept¬ 
ing  estimates  that  are  grossly  in  error,  we  need  tests  that  will 
reject  most  of  the  solutions  that  are  obtained  when  the  target 
maneuvers  radically,  and  hopefully  not  reject  too  many  good  solu¬ 
tions.  Tests  of  this  kind,  using  the  Ekelund  range  estimator  for 
the  last  two  tracking  legs  and  a  three-bearing-rate  range  esti¬ 
mator  (equation  (23)  following)  for  the  last  three  tracking  legs, 
have  been  developed  and  tested  in  reference  (c).  These  tests  will 
be  discussed  briefly  after  the  derivations  for  three  tracking  legs 
are  made. 
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Estimates  From  Three  Tracking  Legs 

Assume  that  three  tracking  legs  are  used  and  we  obtain  two  bear¬ 
ing  estimates,  or  observations,  on  each  leg.  From  these  six 

bearings  we  can  derive  an  exact  range  estimator  that  has  the 

(2)  (3) 

desirable  characteristics  of  the  inexact  estimators  R.  and  R; 

4  4 

in  equations  (2)  and  (3).  That  is,  by  using  six  bearings,  instead 

(2)  (3) 

of  four,  we  can  eliminate  the  biases  in  R.  and  R,  from  formula 

4  4 

inaccuracies,  without  the  necessity  of  using  a  sensitive  estimator 
of  type  R^^  in  equation  (1). 


We  will  use  special  cases  of  three  general  equations,  as  follows: 
rj  Cjk  "  ri  Cik  +  (urk  "  vrkJ  tij  (8) 

°bk  =  vbk  +  rj  V^k  (9) 

“bj  "  Cjk  +  urk  Sjk  (10) 

We  note  that  we  can  replace  the  letter  u  by  the  letter  v  in 
equation  (10) . 


Assume  that  we  wish  to  find  estimators  for  r,  and  u  First,  we 

6  r6 

use  the  second  leg  with  indices  3  and  4  and  the  third  leg  with 
indices  5  and  6.  Special  cases  of  equations  (8),  (9),  and  (10) 
are 


r6  " 

r5  ' 

:56  + 

<“rt 

-  vr6>  C56 

(11) 

r6  " 

r3  ' 

"'36  + 

(ur6 

•  vr6>  C36 

(12) 

“b6  = 

Vb6 

+  r5 

=56' 

'56 

(13) 

V  " 

Vb4 

+  r3 

V 

t34 

(14) 

“b4  " 

"b6 

C.,  +  u  , 
46  r6 

S46 

(15) 

When  we  eliminate  u^,  u^,  r^»  ana  r^  from  these  five  equations 
we  obtain  the  equation 

G,  r  +  H  u  =  M  (16) 

16  1  r6  1 
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where 


G1  =  356  C36  C46/C56  "  8  34’  H1  =  S45  C36/C56  +  8  34  *36  (17) 

(18) 

(19) 


M1  "  Vb4  C36  "  Vb5  C36  C46/C56  +  Vr6  *36  834 


Sjk  "  Sjk7tjk 

If  we  replace  the  subscript  3  by  1  and  4  by  2,  equations  (16), 
(17) ,  and  (18)  become 


G2  r6  +  H2  ur6  =  M2 


(20) 


G2  =  S56  C16  C26/C56  “  812  *  H2  *  S25  C16/C56  +  812  C16  (21) 


M2  “  vb2  C16  '  vb5  C16  C26/C56  +  Vr6  C16  812 


(22) 


where,  now,  6  =  S  /t  . 

We  now  can  solve  equations  (16)  and  (20)  simultaneously  for  r, 

6 

and  ur6-  Using  capital  letters  to  indicate  that  the  bearings 

and  own-ship  velocity  components  are  subject  to  errors,  the 

estimators  for  r,  and  u  ,  are 
t  ro 

R6  "  (M1  H2  ‘  M2  Hi)/(Gl  H2  "  G2  Hl)  (23) 

Ur6  =  (G1  M2  "  G2  Ml)/(Gl  H2  "  G2  V  (24) 

The  estimator  R,  in  equation  (23)  has  formula  accuracy,  but  is 

°  (1) 
much  less  sensitive  to  errors  than  is  the  exact  estimator  R, 

4 

in  equation  (1). 


Now  let  t^  approach  t2»  t^  approach  t^,  and  t^  approach  t^.  Then 
the  G,  H,  M  factors  in  our  estimators  (23)  and  (24)  become 


G 

M 


1 

1 


B,  C.  -  B.,  H.  -  S.,  C,,  +  B.  t.,, 
6  4o  4’  1  46  46  4  46 


Vb4  C46  ”  Vb6  C46  +  Vr6  C46  B4 


(25) 
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G2  “  B6  C26  "  B2*  H2  “  S26  °26  +  B2  *26’ 
M2  ‘  V  C26  -  vb6  C26  +  vr6  C26  *2 


(26) 


where  B~,  B. ,  B.  are  estimates  of  the  bearing  rates  at  times  t  , 
2  4  5  L 

V  V 


Equations  (23) (24) (25) (26)  are  equivalent  to  equations  (A- 30) 
(A-31) (A-32)  of  reference  (c),  where  the  three-bearing-rate  form 
was  derived  directly,  not  via  the  six-bearing  form.  If  it  is 
easier  to  estimate  bearing  rates  and  and  own-ship  velocity  com¬ 
ponents  at  times  t2,,  t^ , ,  t^f  that  are  slightly  less  than  the 

times  t_,  t.  ,  t,  at  the  ends  of  the  tracking  legs,  the  corre- 
L  4  o 

sponding  estimators  (23)  and  (24)  apply  to  time  t&l.  They  will 

have  a  slight  bias  of  the  first  kind  when  used  as  estimators  of 

r,  and  u  ...  If  desired,  the  bias  in  the  range  estimator  could 
o  ro 

be  removed  by  using  an  equation  of  type  (8),  namely, 

R6  "  ^C66'^  ^R6’  +  ^Ur6 '  “  Vr6'^  t6'6^  ('27^ 

but  the  gain  in  formula  accuracy  probably  is  outweighed  by  the 
error  introduced  by  the  use  of  U^,.  The  estimator  R^,  is  recom¬ 
mended,  instead  of  R,  from  equation  (27),  when  it  is  desirable  to 

0 

compute  the  bearing  rate  on  the  last  leg  at  a  time  t^,  other  than 


We  can  use  the  estimators  (23)  and  (24)  to  get  a  complete  set 
of  estimators  for  a  TMA.  Some  TMA  methods  and  the  corresponding 
estimators  are  as  follows: 


TMA  No.  2 


R,  from  equation  (23),  using 
o 

the  G,  H,  M  factors 


the  six-bearing  forms  for 


Ur6'° 

Ub6  =  Vb6  +  R6  S56/t:56 
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TMA  No.  3 


Rg  from  equation  (23),  using  the  three-bearing-rate 
forms  for  the  G,  H,  M  factors 


U  ,  —  U 

r6 

Ub6  “  vb6  +  R6  R6 


TMA  No.  4 


Same  as  TMA  No.  2,  except  that  U  ^  is  obtained  from 
equation  (24),  using  the  six-bearing  forms. 


TMA  No.  5 

Same  as  TMA  No.  3,  except  that  U  ^  is  obtained  from 
equation  (24) ,  using  the  three-bearing-rate  forms 


In  TMA  No.  2  and  TMA  No.  4  a  slightly  better  estimator  of  u 

b6 

is  obtained  by  using  an  equation  of  type  (8)  to  obtain  an 

estimator  for  r  to  replace  R  in  the  equation  for  U,  ,,  since 
jo  bo 

%6  '  vb6  +  r5  V£56 

However,  the  gain  is  slight  and  doesn't  justify  the  additional 
complication. 

TMA  No.  3  and  TMA  No.  5  have  been  studied  in  reference  (c) ,  and 
have  been  compared  with  respect  to  their  effectiveness  in 
corrected-intercept  control  of  the  Mk48  torpedo.  There  is 
little  difference  in  effectiveness,  and  that  is  in  favor  of  TMA 
No.  3  for  the  assumed  error  parameters. 

Discussion  of  Estimators  From  Three  Tracking  Legs 

The  estimators  R,  and  U  ,  in  equations  (23)  and  (24)  are  exact 
o  rb 

and  are  less  sensitive  to  bearing  errors  than  the  corresponding 

12 


Arthur  D  Little  Inc 


exact  estimators  from  two  legs.  The  sensitivity  properties  of 

Rg  are  close  to  those  of  the  inexact  range  estimators  foi  two 

legs.  When  alternate  lag-  and  lead-legs  are  used  for  tracking, 

the  value  of  |  Hj  |  usually  is  an  order  of  magnitude  greater 

than  j  H^| ,  thereby  eliminating  a  major  source  of  sensitivity. 

Some  results  for  R,  with  the  three-bearing-rate  forms  of  G,  H,  M 
0 

are  given  in  reference  (c) . 

The  problem  of  estimating  u^  from  bearing  observations  alone  is 
a  difficult  one.  After  a  careful  study  of  the  problem  we  have 
come  to  the  conclusion  that  the  estimator 

Ur6  *  0  (28) 

is  hard  to  beat.  By  accepting  a  100  percent  bias  we  avoid  all 
other  biases,  some  of  which  can  easily  be  greater  than  100  per¬ 
cent.  We  also  have  zero  variance.  The  estimator  (24)  is  quite 
sensitive  to  bearing  errors  and  to  non-linear  target  motion, 
although  it  is  less  sensitive  to  these  variations  than  other 
estimators,  such  as  CHURN,  that  we  have  studied. 

In  some  applications  the  range  component  of  target  velocity  is 
not  needed,  or  it  plays  a  minor  role.  For  example,  the  lead 
angle  for  interception  of  the  target,  assumed  to  be  in  linear 
motion,  doesn't  require  the  range  component  of  target  velocity, 
if  the  guide  point  is  on  the  current  bearing  line.  Also,  any 
realizable  value  of  the  range  component  of  target  velocity  has 
little  effect  on  the  lead  angle  when  the  guide  point  is  close  to 
the  current  bearing  line.  Hence,  under  these  conditions  any 
estimator  that  doesn't  yield  impossible  values  can  be  used.  And 
impossible  values  can  be  avoided  when  using  the  estimator  (24)  by 
a  simple  acceptance/rejection  test  based  on  the  magnitude  of  the 
estimate.  A  simpler  way  to  avoid  impossible  values  is  to  use  the 
estimator  (28),  which  proved  to  be  slightly  superior  to  the  esti¬ 
mator  (24)  with  an  acceptance/rejection  test  in  some  fire-control 
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problems,  as  reported  in  reference  (c).  This  conclusion  does 
not  necessarily  apply  to  otner  uses  of  the  u^  estimators. 

Acceptance/Rejection  Tests 

The  use  of  three  tracking  legs  requires  a  tracking  interval 
approximately  50  percent  larger  than  that  required  for  two  track¬ 
ing  legs,  and  provides  some  advantages.  First,  it  enables  us  to 
obtain  a  comparatively  rugged  range  estimator  (23)  that  has 
formula  accuracy.  Second,  it  provides  an  estimator  (24)  of  the 
range  component  of  velocity  that  might  be  competitive  with  the 
simple  estimator  (28)  in  some  applications.  Third,  it  provides 
us  with  a  means  of  constructing  simple  acceptance/rejection  tests 
to  avoid  gross  errors  produced  by  non-linear  target  motion.  In 
our  opinion,  it  is  only  the  third  advantage  that  justifies  the 
use  of  more  than  two  tracking  legs. 


Some  acceptance/rejection  tests  are  described  in  reference  (c) , 
together  with  the  results  obtained  from  a  computer  simulation. 
Three  tracking  legs  are  run  and  the  test  is  made.  If  the  out¬ 
come  is  positive,  the  solution  is  accepted.  If  the  outcome  is 
negative,  the  solution  is  rejected  and  another  tracking  leg  is 
run,  continuing  until  the  first  positive  outcome.  A  simple  test 
of  this  type  for  the  range  estimate  is  the  following: 

Let 

R  =  Ekelund  range  estimate  from  the  last  two 
E 

tracking  legs,  using  an  estimator  of  type  (6) 

in  which  the  v,  and  B  components  are  esti- 
b 

mated  for  the  end  of  the  tracking  legs 

Rd  =  range  estimate  from  the  last  three  tracking 
legs,  using  an  estimator  of  type  (23)  with 
the  G,  H,  M  factors  computed  from  the  bearing 
rates  at  the  ends  of  the  last  three  legs 
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. . 


(rrain  »  rmax)  “  acceptable  range  interval 


a  =  coincidence  factor,  0  <  a  <  1 
r  r 


The  test  conditions  are: 


a 

r 


<  Rj/S,,  <  l/ar 


(29) 


r  <  R  < 
min  E 


r 

max 


min 


<  R. 


r 

max 


(30) 


If  all  conditions  are  satisfied,  accept  the  solution;  otherwise, 
reject  it.  Then  there  remains  the  question  of  which  range  esti¬ 
mate,  R  or  R  .  to  use  when  the  solution  is  accepted.  This 
question  was  made  a  part  of  the  simulation  study  in  reference 

(c) ,  as  was  the  question  of  what  values  to  use  for  a  ,  r  , 

r  min 

and  r 

max 


At  first,  only  condition  (29)  was  used.  Conditions  (30)  were 
added  to  eliminate  a  few  very  wrong  solutions  for  which  condi¬ 
tion  (29)  was  satisfied.  For  this  purpose,  it  was  found  that 

the  values  for  r  .  and  r  could  be  chosen  to  be  any  reason- 
min  max 

able  limits  within  which  the  range  could  lie  and  the  desired 

(fire-control)  objectives  could  be  achieved.  The  bias  in  R„ 

E 

from  formula  inaccuracy  made  it  necessary  to  use  what  might 
appear  to  be  small  values  of  a^,  to  avoid  rejecting  too  many 
"good"  solutions. 


Other  acceptance/rejection  tests  were  tried,  including  one  that 
seldom  accepts  a  poor  solution.  However,  it  requires  four 
tracking  legs  and  often  rejects  good  solutions.  A  simple  test 
of  the  form 


U  <  u 
r'  max 


(31) 


was  used  when  the  estimator  (24)  was  used  for  the  range  component 
of  velocity. 
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Time  Correction 


This  method  of  range  estimation  is  described  in  reference  (b) . 

A  critique  of  the  method  was  made  in  reference  (a),  and  a  reply 
was  made  in  reference  (d) .  We  believe  that  a  fair  and  impartial 
summary  of  the  method  is  as  follows: 

•ff 

Time  correction  is  a  method  of  obtaining  an  estimator  R  for  the 
* 

range  at  a  time  t  that  is  computed  from  the  bearing  observa¬ 
tions  and  the  corresponding  times  at  which  the  observations  are 

*  — 
made.  The  equation  for  R  contains  two  additional  bearings  B 
* 

and  B  that  must  be  estimated  by  interpolation  or  extrapolation 

from  the  bearing-time  plot;  their  use  can  be  avoided  by  using 

small-angle  approximations,  the  validity  of  which  can't  be 

* 

stated  without  estimating  the  two  bearings.  The  time  t  at 

* 

which  the  estimator  R  applies  is  only  partially  under  the 
control  of  the  tracking  submarine;  it  is  unlikely  to  be  close 
to  the  time  at  which  a  range  estimate  is  needed. 

In  reply  to  the  question  of  how  one  would  get  an  estimate  of 

* 

the  range  at  a  time  different  from  t  ,  it  is  stated  in  reference 
(d)  that  a  second  time-corrected  range  estimate  would  be  made 
and  used  to  get  a  complete  TMA,  from  which  the  desired  range 
estimate  presumably  would  be  obtained.  And  we  are  referred  to 
Chapter  5  of  reference  (b)  for  the  details. 

A  method  of  getting  a  complete  TMA  from  four  bearings,  taken  on 
three  tracking  legs,  is  described  in  Sections  5. A  and  5.5  of 
reference  (b).  In  our  notation,  the  epoch  t^  occurs  on  the 
first  leg,  t2  and  t^  on  the  second  leg,  and  t^  on  the  third  leg, 
according  to  Figure  5-4.  In  our  opinion,  this  method  has  sev¬ 
eral  undesirable  properties.  Since  t^  and  t^  occur  on  different 
tracking  legs,  S^  will  have  a  large  bias  from  delta  biases  in 
the  bearings;  similarly,  for  S^.  These  biases  can  be  reduced 
significantly  by  pairing  bearings  on  the  same  tracking  leg,  as 
explained  earlier.  Secondly,  the  four-bearing  TMA  obtained  by 
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the  method  described  in  reference  (b)  is  equivalent  to  (that  is, 
identical  to,  in  the  absence  of  errors)  the  ordinary  four-bearing 
TMA,  as  stated  on  page  5.8  of  reference  (b) .  Hence,  it  is  very 
sensitive  to  random  bearing  errors  and  to  small  deviations  from 
linear  target  motion.  If  we  are  willing  to  accept  these  sensi¬ 
tivities,  why  not  use  the  simpler,  and  equivalent,  direct  solu¬ 
tion,  for  which  in  equation  (1)  is  the  range  estimator? 

We  believe  that  a  better  method,  if  it  is  desired  to  use  time 
correction,  is  the  following:  Take  two  bearing  observations 
on  each  tracking  leg  and  pair  them  to  avoid  large  biases  from 

•k 

delta  biases.  Use  legs  2  and  3  to  get  the  estimator  R1  at  time 

it  it  * 

t^,  and  legs  1  and  3  to  get  the  estimator  at  time  t^.  To 

get  estimators  for  r,  and  u  , ,  for  example,  use  an  equation  of 

type  (8)  with  t.  ■  t- ,  t  =  t0,  t,  «  t,  and  solve  for  u  , . 

i  l  j  2  k  6  *  r6 

Then  use  another  equation  of  type  (8)  with  t^  =  t^,  t  =  t^  =  t^ 
and  solve  for  r^. 

We  submit  that  the  above  TMA  method,  in  the  absence  of  errors, 
is  identical  with  the  direct  method  that  resulted  in  the  esti¬ 
mators  R,  and  U  ,  in  equations  (23)  and  (24).  The  time-correction 
o  ro 

method  has  all  the  errors  of  the  direct  method,  plus  a  few  more. 

The  additional  errors  are  those  that  occur  in  estimating  the  four 

_  £  _  £  _  £  _ 

bearings  B^,  B^,  B^,  and  B2  at  the  computed  times  t^,  t^,  t2»  and 

t2  from  the  bearing-time  plot  by  interpolation  or  extrapolation; 

or  alternatively,  the  biases  that  occur  from  the  use  of  the  small- 

angle  approximations.  Is  there  any  gain  that  offsets  these  errors 

and  the  additional  complications?  We  believe  that  the  direct 

method  is  preferable,  and  that  the  very  simple  estimator  (28)  that 

is  used  in  TMA  No.  2  and  TMA  No.  3  may  be  even  better  for  some 

purposes  than  the  estimator  (24). 

The  CHURN  TMA 

The  CHURN  TMA  uses  estimators  for  the  rectangular  coordinates 
(xq,  yQ)  of  the  relative  position  of  the  target  at  the  start  of 
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tracking  and  for  the  rectangular  components  (u  ,  u  )  of  target 

x  y 

motion.  The  estimators  are  obtained  by  minimizing  a  sum  of 
(perhaps  weighted)  squares  of  residuals  measured  normal  to  the 
bearing  lines. 

The  results  of  a  study  of  the  biases  in  the  CHURN  TMA  are  given 
and  discussed  in  reference  (c).  Delta  biases  in  the  bearing 
errors  produce  large  errors  in  the  values  of  the  position  and 
motion  parameters.  Also,  random  bearing  errors  produce  very 
large  biases  in  the  parameters  under  some  conditions,  which  was 
not  expected.  Large  deviations,  and  even  moderate  deviations, 
from  a  linear  course  produce  large  biases  in  the  parameters. 

The  stochastic  bias  and  variance  were  studied  by  Monte  Carlo 
simulation.  The  very  large  biases  from  strictly  linear  motion 
and  random  bearing  errors  were  so  unexpected  that  they  were 
ascribed  to  undetected  errors  in  the  computational  procedure 
until  they  were  understood  and  could  be  explained.  When  the 
mechanism  that  produces  the  biases  was  understood,  a  simple 
algorithm  could  be  constructed  to  compute  the  stochastic  bias 
without  recourse  to  Monte  Carlo  simulation. 

The  stochastic  bias  occurs  because  a  radial-coordinate  residual 
is  used  in  a  rectangular  coordinate  system,  thereby  producing 
estimators  for  the  parameters  that  are  far  from  linear  in  the 
bearings.  The  sum  of  squares  of  the  residuals  can  be  written 
as  a  sum  of  two  terms,  one  arising  from  biases  in  the  solution 
and  the  other  from  the  random  variations  in  the  bearings.  The 
minimization  of  the  sum  of  the  two  terms  produces  a  biased 
solution  that  has  a  much  smaller  value  than  the  unbiased  solu¬ 
tion  for  the  random  component  of  the  sum  and  a  non-zero  value 
for  the  bias  component  of  the  sum.  The  bias  in  the  solution  is 
in  the  direction  of  predicting  that  the  target  track  over  the 
tracking  interval  lies  closer  to  the  tracking  submarine  than 
does  the  actual  target  track. 
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The  magnitude  of  the  stochastic  bias  depends  on  several  factors, 
of  which  the  duration  of  the  tracking  interval  and  the  standard 
deviation  in  the  bearing  errors  appear  to  be  the  most  important. 
For  tracking  intervals  of  10  minutes  or  less  and  standard  devia¬ 
tions  of  0.5  degrees  or  more,  the  stochastic  bias  is  so  large  as 
to  make  the  solution  of  doubtful  value. 

Large  deviations  from  a  linear  course,  initiated  near  the  middle 
of  the  tracking  interval,  produce  large  errors  in  the  parameters. 
In  general,  CHURN  is  slow  to  respond  to  changes  in  target  course 
or  speed  after  tracking  has  been  in  progress  for  a  long  period. 
Thus,  changes  in  course  or  speed  that  occur  late  in  the  tracking 
interval  will  have  little  effect  on  the  solution.  This  charac¬ 
teristic  of  CHURN  is  an  important  one  in  fire-control  applica¬ 
tions. 

A  change  in  target  course  or  speed  late  in  the  tracking  interval 
often  produces  a  large  bias  in  the  lead  angle.  The  estimate  of 
range  may  not  be  greatly  in  error,  but  the  estimate  of  the  veloc¬ 
ity  component  u^  normal  to  the  bearing  line,  from  which  the  lead 
angle  is  computed,  will  not  reflect  the  current  target  motion. 

The  large  lag  in  the  lead  angle  from  CHURN  is  intolerable  in 
post-launch  control  of  a  wire-guided  torpedo  in  the  corrected- 
intercept  mode.  The  sole  purpose  of  post-launch  tracking  and 
control  is  the  detection  of  changes  in  target  motion  that  occur 
after  launch,  and  in  time  to  permit  corrections  to  the  torpedo 
course  to  be  computed  and  applied.  For  this  purpose  we  believe 
that  TMA  No.  3  and  T&A  No.  5 — and,  we  prefer  TMA  No.  3 — are 
superior  to  CHURN  by  a  wide  margin.  If  a  moderately  good  esti¬ 
mate  of  the  bearing  rate  for  the  current  target  motion  can  be 
obtained  from  the  bearing-time  plot  or  by  the  exponential¬ 
smoothing  method  used  in  reference  (c) ,  the  estimate  of  the 
normal  velocity  component  will  be  accurate  enough  to  yield  a 
good  value  for  the  lead  angle. 
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Filtering  and  Smoothing  Techniques 


The  purpose  of  filtering  and  pre-smoothing  techniques  is  to  re¬ 
move  the  random  bearing  errors,  to  the  extent  that  this  is  pos¬ 
sible.  For  the  sake  of  argument,  assume  that  the  random  bearing 
errors  can  be  reduced  by  several  orders  of  magnitude  through  the 
use  of  a  sophisticated  technique,  such  as  Kalman  Filtering,  and 
effectively  eliminated  thereby.  What  effects  will  this  change 
have  on  the  total  error  structure  of  the  estimators? 


Only  the  stochastic  bias  and  the  variance  are  affected.  The 
biases  from  formula  inaccuracies,  from  biases  in  the  bearing 
errors,  and  from  non-linear  target  motion  are  unaffected.  The 
net  effect  on  the  total  bias  may  be  large  or  insignificant, 
depending  on  the  magnitudes  of  the  other  biases  relative  to  that 
of  the  stochastic  bias.  Only  when  the  stochastic  bias  is  the 
dominant  bias  will  the  net  effect  on  the  total  bias  be  large. 
This  dominance  will  occur  in  some  cases  with  the  estimators 
obtained  from  CHURN,  particularly  when  the  target  motion  is 
nearly  linear  and  the  tracking  interval  is  short.  We  believe 
that  the  CHURN  TMA  is  a  likely  candidate  for  the  application  of 
filtering  techniques. 


A  larje  reduction  in  the  random  bearing  errors  may  reduce  the 
variances  in  the  estimators  from  this  source  to  the  point  that 
variations  in  other  estimates,  such  as  own-ship  velocity  compo¬ 
nents,  are  dominant.  The  net  effect  of  a  large  reduction  in 
the  variance  of  an  estimator  depends  on  the  application  and  the 
magnitude  of  the  total  bias.  In  fire-control  applications  the 
optimal  variance  usually  is  not  zero  when  the  total  bias  is 
large  in  magnitude,  in  which  case  a  large  reduction  in  the 
estimator  variance  may  produce  a  net  loss. 

Another  possible  gain  from  the  use  of  filtering  techniques  is 
that  of  obtaining  more  accurate  estimates  of  bearing  rates  than 
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can  be  obtained  from  a  smoothed  bearing-time  plot  or  by  exponen¬ 
tial  smoothing,  for  use  in  our  TMA  No.  3  or  TMA  No.  5.  The  need 
for  accurate  estimates  of  bearing  rate  is  particularly  great  in 
the  corrected-intercept  control  mode  that  we  proposed  in  refer¬ 
ence  (e)  for  use  with  degraded  tracking.  This  use  of  filtering 
techniques  may  prove  to  be  more  important  than  other  uses  that 
have  been  proposed. 

Test  and  Comparison  Program 

Some  methods  of  computing  biases  and  variances  of  the  estimators 
are  described  in  reference  (c).  Ue  propose  that  these  methods 
be  reviewed,  revised  as  necessary,  and  applied  to  the  principal 
estimators  and  tests  described  above  and  to  others  that  appear 
to  be  competitive. 

Computer  simulation  has  been  found  to  be  the  simplest  method  of 
computing  biases  from  non-linear  target  motion.  This  procedure 
also  can  be  used  to  compute  estimator  biases  from  formula  inac¬ 
curacies  and  from  bearing  biases,  although  simple  analysis  also 
can  be  used.  The  point  here  is  that  the  use  of  computer  simula¬ 
tion  for  the  computation  of  biases  from  all  three  sources  can  be 
done  with  little  more  effort  than  that  needed  for  the  non-linear 
motion. 

Biases  and  variances  from  random  bearing  errors  can  be  handled 
by  simple  analysis  in  many  cases.  We  had  to  resort  to  Monte 
Carlo  simulation  in  the  use  of  CHURN,  and  it  may  be  necessary 
for  some  other  estimators. 

The  analysis  of  errors  should  be  extended  to  include  errors  in 
other  parameters,  such  as  own-ship  position  and  velocity  compo¬ 
nents.  The  analysis  also  should  include  tests  of  the  accept/reject 
type. 
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